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Abstract 

We consider dynamical systems associated to Lax pairs depending rationnally on 
a spectral parameter. We show that we can express the symplectic form in terms of 
algebro-geometric data provided that the symplectic structure on L(X) is of Kirillov 
type. In particular, in this case the dynamical system is integrable. 
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1 Introduction. 



The algebro-geometrical approaches to integrable systems are the most powerful ways 
of solving the Lax equations of motion. These methods allow to introduce natural coor- 
dinates on phase space, namely g action variables which characterize the spectral curve 
of genus g and a divisor of g points on this curve which is equivalent to a point on the 
Jacobian and can be related to angle variables. These coordinates are also known as the 
separated Sklyanin's variables ||. It is an important and natural problem to express the 
symplectic form in terms of these coordinates. The first conjecture was made by Veselov 
and Novikov in Jl|. More recently the question has been reconsidered by Krichever and 
Phong [|J who proposed a general method to prove such a connection. Application of 
their technique to the spin Calogero model @, [|, U met some difficulties and led to im- 
provements and simplifications of the method. In particular the role of the branch points 
has been recognized in ||. The spin Calogero model however is rather specific and we 
propose here to analyze the situation corresponding to a generic rational Lax pair. In 
this case the natural symplectic structure is a coadjoint orbit one, and we find that we 
have to use it in all details in order to complete the calculation of the symplectic form. 
This shows the nice interplay between the group theoretical and analytic approaches to 
integrable systems. 

2 The setup. 

Let us consider a N x N Lax matrix L(X), depending rationally on a spectral parameter 
A 6 C, with poles at points Xk 

L(\) = L + J2 L k (1) 

k 

where L = Diag(ai, • • • , a^v) is independent of A and L k is the polar part of L(X) at Xk, 
ie. L k = J2r=-n k ^fc,r(A — Afc) r . We assume that L k lives in a coadjoint orbit of the group 
of N x N matrix regular in the vicinity of A = A^, i.e. 

L k = (gkA^ 1 ^ 

where Ak(X) is a diagonal matrix with a pole of order at A = Xk, and has a regular 
expansion at A = Xk- The notation ()_ means taking the singular part at A = Xk- This 
singular part only depends on the singular part (Ak)- and the first rik coefficients of the 
expansion of gu in powers of (A — Xk). We see that (Ak)- specifies the coadjoint orbit, 
and that the physical degrees of freedom are contained in the first n& coefficients of gk(X). 
Since Ak commutes with diagonal matrices one has to take the quotient by — > g^dk 
where (4(A) is a regular diagonal matrix, in order to correctly describe the dynamical 
variables on the orbit. 

To interpret L as living in a coadjoint orbit we consider the direct product of loop 
groups defined for each A&. Its Lie algebra is a direct sum of Lie algebras Qk where an 
element of Qk is a regular sum X k (X) = Z) n >o A fcn (A — Xk) n and the Lie bracket is such 
that [Xk(X), Xi(\)] = if k ^ I. The dual of this Lie algebra can be viewed as iV x JV 
matrix valued meromorphic functions L(\) with poles at the A& and vanishing at oo. The 
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duality is expressed by: 

< L(X),® k X k (X) >= ]TRes Afc Tr(L(A)X fc (A)) 

k 

Note that L can be uniquely decomposed as J2k ^k, where L k is the polar part at and 
that only L k appears in the residue ResA fc . One can add an invariant character Lq to L. 

The dimension of the orbit of L k is N(N—l)n k so that L(X) depends on J2k N(N—l)nk 
degrees of freedom. On this phase space we have a standard Kirillov symplectic form: 

u = Res Afc Tr((A k )_g^5g k A g k l 5g k ) d\ (2) 

k 

In the algebro-geometric approach we study the eigenvector equation: 

(L(A) - Ail) *(A, a*) = (3) 

and view the eigenvector \I/(A, fi) with eigenvalue \i as an analytic section of a natural line 
bundle on the Riemann surface defined by the characteristic equation: 

T : R(X, n) = det(L(A) - fi 1) = (4) 

If N is the dimension of the Lax matrix, the equation of the curve is of the form: 

AT-l 

T : J2(A,Ai) = (-0)*+£ r «( A K = O (5) 

9=0 

The coefficients r q (X) are polynomials in the matrix elements of L(X) and therefore 
have poles at A&. Since the Lax equation L = [L, M\ is isospectral, these coefficients are 
time-independent and are related to the action variables. 

From eq.(||D, we see that the spectral curve appears as an iV-sheeted covering of the 
Riemann sphere. To a given point A on the Riemann sphere there correspond N points 
on the curve whose coordinates are (A,/ii), ■ • • (A, a*at) where the /ij are the solutions of 
the algebraic equation i?(A,/j) = 0. By definition are the eigenvalues of L(X). 

We assume for simplicity that all the a,'s are different. Then on the spectral curve, we 
have iV points Qi = (A = oo,/Xj = Oj) above A = oo. The analyticity properties of L{X) 
are invariant under conjugation by constant matrices. To preserve the normalization at 
oo these matrices have to be diagonal. Generically, these transformations form a group 
of dimension N — 1 and we will have to factor it out. 

Before doing complex analysis on T, one has to determine its genus. A general strategy 
is as follows. As we have seen, T is a iV-sheeted covering of the Riemann sphere. There is 
a general formula expressing the genus g of an iV-sheeted covering of a Riemann surface 
of genus go (in our case go = 0). It is the Riemann-Hurwitz formula: 

2g - 2 = N(2g -2) + v (6) 

where v is the branching index of the covering. The branch points occur at the zeroes of 
d^R. The number of its zeroes is the same as the number of its poles, which are located 
above the points X k . One gets: 

N(N-l) ^ 
1 k 
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For completeness of the method it is important to observe that the genus is equal to 
the number of action variables occurring as independent parameters in the eq.fl^) which 
should also be half the dimension of phase space. This phase space Ai is the above 
coadjoint orbit, quotiented by the action of constant diagonal matrices. The orbits of this 
action are of dimension (JV — 1), since the identity does not act. One has to perform a 
Hamiltonian reduction by this action. First one fixes the momentum, yielding (JV — 1) 
conditions, and then one takes the quotient by the stabilizer of the momentum which 
is here the whole group since it is Abelian. Hence the dimension of the phase space is 
reduced by 2(JV — 1), yielding: 

dimM = (JV 2 - N)J2n k -2(N - 1) = 2g 

k 

Let us now count the number of independent coefficients in eq.(|5]). It is clear that rj(X) 
is a rational function of A. The value of rj at oo is known since fij — > a,. Moreover Tj(X) 
has a pole of order jn k at A = X k . Hence it can be expressed on j J2 k n k parameters namely 
the coefficients of all these poles. So we have altogether ^N(N + 1) J2k n k parameters. 
They are not all independent however. Indeed above X = X k the various branches can be 
written: 

where all the coefficients of , • • • , are fixed and non-dynamical because they are the 
matrix elements of the diagonal matrices (-A&)_, while the regular part is dynamical. 
This implies Nn k constraints on the coefficients of r 3 - (note that rj is the symmetrical 
function (Jj(/J,i, ■ ■ ■ , hence the n k highest order terms in Tj are fixed). We are left with 
|iV(iV — 1) J2k n k parameters, that is g + iV — 1 parameters. 

It remains to take the quotient by the action of constant diagonal matrices. Consider 
the Hamiltonians H n = (1/n) resA =00 Tr (L n (A)) dX, i.e. the term in 1/A in Tr(L n (A)). 
These are functions of the Tj(X). We show that they are the generators of the diagonal 
action. First we have: 

ResA=ooTr (L n (X))dX = n Res A=00 Tr {L^ 1 ]T L k {X))dX 

k 

= nRes A =ooTr(L^ 1 L(A))rfA 

since all L k are of order 1/A at oo. Computing with the above Kirillov bracket one obtains 
{H n , L(n)} = — [L'q^ 1 , L{jj)\ which is the coadjoint action of a diagonal matrix on L(fi). 
Since L is generic the Lq generate the space of all diagonal matrices, so we get exactly 
N — 1 generators Hi, ■ ■ ■ , i^jv-i- 111 the Hamiltonian reduction procedure, the H n are the 
moments of the group action and are to be set to fixed (non-dynamical) values. Hence 
when the system is properly reduced we are left with exactly g action variables. 

The eigenvector line-bundle has Chern class — (g + N — 1). We get a non- vanishing 
section by requiring that the eigenvector ^f(P) has its first component equal to 1. Then 
it has (g + N — 1) poles on the spectral curve 0. Since L is diagonal at oo the natural 
eigenvectors of L(X) above A = oo are the canonical basis vectors. Since we impose that 
^1 = 1 this implies that ^(P) has poles at (JV — 1) points above A = oo, hence we are 
left with g dynamical poles. Note that this is the number of angle variables and we are 
led to describe the system in terms of the g action variables entering the equation of the 
spectral curve and g points on this curve, namely the dynamical poles of 
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3 Symplectic form. 



We have seen that the Lax pair description of a dynamical system naturally provides 
coordinates on phase space, namely g independent action variables Fi which parameterize 
the spectral curve F, and g points ^ = (A,, /ij) on the spectral curve, which we called 
the dynamical divisor. It is important to express the symplectic form in terms of these 
coordinates. The phase space appears as a fibered space whose base is the space of moduli 
of the spectral curve, explicitly described as coefficients of the equation R(X, u) — of 
the spectral curve, and the fiber at a given R is the Jacobian of the curve R = 0. On 
this space we introduce a differential 5 which varies the dynamical variables Fj, Aj, /ij 
subjected to the constraint R{ Fi y(Xi, fJ>i) = 0. 

We will need an auxiliary fiber bundle above the same base whose fiber is T x Jac(r). 
We extend S to this space by keeping the previous definition on the Jac(T) part and on the 
r part, we differentiate any function of Fi, A, /i with Rsp.\(X, u) = keeping A constant. 
Remark that A is universally defined on the whole family of curves. For a function /(P; Pj) 
if we take A as a local parameter 5f = J2i d F J8Fi. At a branch point the local parameter 
is \x and we have: 

5f = d,f5fi + J2d F J5F l , with 6fM=- 1 £ d Fi R {Fi} (\, y)8Fi (7) 

Note that at a branch point d IJi RjF.\(X, u) = 0, hence the differential Sf acquires a pole 
even though / is regular. Remark however that if / depends rationally on A and the Pj, 
5f is regular at the branch points. 

At each point P(A, a) on T is defined a column eigenvector ^(P) of the Lax matrix up 
to normalization. This allows to define a matrix \&(A) whose columns are the N vectors 
ty(Pi) at the iV points Pi above A (assuming we have chosen locally some ordering of the 
sheets). We also consider its inverse matrix \1/ _1 (A) and denote its iV lines by ^"^(Pj) 
(for the same ordering of the Pi). In particular we have < \I/( -1 )(P), ^(P) >= 1. Note 
that ^^(P) has poles at the branching points of the covering (A,/i) — > A since the 
determinant of ^ vanishes here. 

We define a three-form K on our extended fiber bundle, and we regard it as a one-form 
on T whose coefficients are two-forms on phase space. 

K = K 1 + K 2 + K 3 (8) 

K x = < ¥- 1) {P)SL(X) A5^(P) > dX 

K 2 = <¥- 1 \P)5/iA5^(P)>dX 

K 3 = 5(log<9 M P) A5/jdX 

Proposition. Let us define the two-form on phase space: uo = J2k,i R- es P fc t K where P k>i 
are the points above the poles A& of L(X). Then we have: 

9 

uj = 2 5Xi A 5 Hi 

i=l 

hence uj is a symplectic form on phase space. 

Proof. The sum of the residues of K seen as a form on F vanishes. The poles of 
K are located at four different places, first the dynamical poles of then the poles at 
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the Pfc,i coming from L and /i, next the poles above A = oo coming from 9 and dX, and 
finally the poles at the branch points of the covering coming from the poles of 9^~^ and 
from eq.(^). 

Let us compute the residues at the dynamical poles (ui, ■ ■ ■ , v g ). We write the coordi- 
nates of these points as: Vi = (Aj, /Xj) for % — 1, • • • g. Near such a point we can choose A 
as a universal local parameter and 9 = 1/(A — Aj) x 9 reg hence: 

69 = ^-^ + 0(1)) 

Since (L - n)9 = and - /j) = 0, we have (5L - 5fi)9 + (L - /i)59 = 0. 

Multiplying by 9 ( -~ 1 ' ) we get 9^^5L9 = Sfi, therefore: 

Res^iTi = 8n\ Vi A 8Xi 

Here 5/j, is to be seen as a meromorphic function on T given by eq.(^). However varying 
R(Xi, fii) = we obtain: 

8/i\ v . = 5/J,i+ ' C A 



5A 3 ; 



and the second term does not contribute to the wedge product. 

The contribution of K 2 is exactly the same. Finally, K 3 is regular at and does not 
contribute to the residue at this point. So we finally get: 

Res Ui K = 25 fa A 5X { (9) 

We now show that there are no residues at the branch points due to the proper choice 
of K2 and K3. Let us look at the term K\. At a branch point b, \E^ _1 ) has a simple pole, 
5L is regular, 59 has a simple pole due to eq.(^) and the form dX has a simple zero, hence 
the considered expression has a simple pole at b. To compute its residue it is enough 
to keep the polar part in 59, i.e. to replace 59 by d^95fi (recall that fj, is a good local 
parameter around b). We get: 

ResfeKi = Res b <9 { ~ 1) 5Ld fl 9 > ASfidX 

= Rest < 9 ( ~ 1) (5L - 5fi)d^9 > ASfidX 

where in the last equation we have used the antisymmetry of the wedge product to replace 
5L by SL — Sfi. Using again the eigenvector equation (L — fi)9 = 0, and varying the point 
(A, /i) on the curve around b one gets 

where dX/d/i vanishes at the branch point. We then differentiate with 5 and get: 
Res b < 9 { ~ 1] (5L - 5fi)d^9 > ASfidX = Res b < 9 { ~ 1] 59 > ASfidX 

-Res b f^vlA > A5fi dX 

\dfi dX I 
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The first term exactly cancels the term Res^-K^- The second term gives a non- vanishing 
contribution 

Res fe — A SpL dX 

fi — fib 

Indeed let us call ( = (dX/du)(dL/dX)ty which vanishes at b = (A&, Then ( = 
(/i — /ib)Ci hence 5( = —5fib/(n — Hb)( + <fy* C2 + C3 with £ 3 regular. The second term does 
not contribute due to the antisymmetry of the wedge product and the third term has no 
residue. Using eq.flTUD we have < ^ _1 ^^x^ >— 1 yielding the above formula. This 
contribution is exactly canceled by the contribution of K3. 

We now compute the residues above A = 00. Recall that we consider a reduced 
Hamiltonian system under the action of diagonal matrices. To fix this action on can 
normalize the eigenvectors at 00 so that ipi{Qj) = X5ij + 0(1) for i, j — 2, ■ ■ ■ , N. Notice 
that L = Lq + 0(1/A) where L$ is non-dynamical so 5Lq = 0, and that /i = + 0(1/ X) 
around Qi hence 5L = Su. — 0(1/A). Moreover vanishes at Qi and dX has a 

double pole. Altogether K\ and K 2 are regular at Qi since (5*f>)(Qi) = 0(1) due to the 
normalization condition. Finally K% is also regular since on the sheet fi = /ij(A) one can 
write d^R = Yij^A^i — Hj) yielding 5\ogd^R = 0(1/A). All this shows that K has no 
residues above A = 00. 

Proposition. The symplectic form u is given by: 

w = 2 £ to, A fy, = 2 £ Res Afc Tr((A fe )-#fc x <^ A g^5g k ) dX (11) 

i=l k 

where (A i; m), % — 1, • • ■ g, are the coordinates of the points of the dynamical divisor D. 

Proof. Let us compute the residues at the poles X k of K\, where only L k contributes. 
Since locally we have L = gkAkg^ 1 around A = X k , we may identify the matrix ^(A) with 
g k . More precisely we have ^(A) = g k d k and \& _1 (A) = d^g^ 1 with d k a diagonal matrix. 
The residues are obtained by integrating over small circles surrounding each of the N 
points P k ,i above Afc. We can choose these small circles so that they project on the base 
A on a single small circle surrounding X k . Then we get 
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EResp,,^ = f <¥- 1 \P l )SL(X)AS^(P l )> dX 

i = l i=l Lm J °k,i 

= — l Tr(^- 1 (A)5L(A) A 5§(A)dA) (12) 
2^7^ Jc k 

where we used the fact that \& _1 (A) is equal to the matrix whose lines are the vectors 
vJm -1 )^). xhe trace has been reconstructed in eq. (jl2"l) because ^(Pi), i = 1, ■ ■ ■ N, form 
a basis of eigenvectors. Using the identification of ^(A) in terms of g k gives: 

Res Afe ^i = Res Afc Tr^ 1 ^ 1 (5g k (A k )_g k 1 - g^A^^g^dg^ 1 ) A 

(5g k d k + g k 5d k ))dX 
= -2 Res Afc Tr((A fe )_^ 1 55f fe A g k l 5g k ) dX 

+Res Afc Tr{g^5g k [(A k )_, Sd k d^}) dX (13) 

The last term vanishes because it involves the commutator of two diagonal matrices. To 
compute the residues of K 2 at X k we remark that 5/i is regular on all sheets above X k . 
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This is because due to the form of L(X), we have fx = (Ak)- + regular. Since (Ak)- 
characterizes the coadjoint orbit and is not dynamical, one has to take S(Ak)_ = 0. It 
follows immediately that K 2 is regular. To compute the residue of K 3 we note that if d^R 
has a pole of some order m at P^i and can be written d^R = c(A) /(A — \k) m where c(A) is 
regular and non vanishing, we get 5 (log d^R) = <51ogc(A) which is regular, since 5\ = 
and SXk = 0. Hence K 3 has no residue again because (5/x is regular. 

This proposition means that the coordinates (Aj,//j) of the point z/j of the dynamical 
divisor are canonical coordinates. This type of result can be obtained in the r-matrix 
approach for specific models like the Toda chain B. 

4 Conclusion 

We have shown under quite general but necessary assumptions that the natural symplectic 
structure on Lax pairs can be expressed in terms of algebro-geometric data. This result 
shows the nice interplay between the analytical and the group-theoretical approaches to 
integrable systems. We are able to show that (Ai,//j) are canonical coordinates only using 
the fact that L parametrizes a coadjoint orbit, specified by constant matrices (Ak)- and 
La. 
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